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Abstract. We report circumstantial evidence for the first de- 
tection of pulsed high-energy 7-ray emission from a millisec- 
ond pulsar, PSR J0218+4232, using data collected with the En- 
ergetic Gamma Ray Experiment (EGRET) on board the Comp- 
ton Gamma Ray Observatory (CGRO). The EGRET source 
3EG J0222+4253 is shown to be spatially consistent with PSR 
J0218+4232 for the energy range 100 - 300 MeV. Above 1 
GeV the nearby BL Lac 3C 66A is the evident counterpart, 
and between 300 MeV and 1 GeV both sources contribute to 
the 7-ray excess. Folding the 100-1000 MeV photons with an 
accurate radio ephemeris of PSR J0218+4232 yields a double 
peaked pulse profile with a ~ 3.5er modulation significance 
and with a peak separation of ~ 0.45 similar to the 0.1-10 
keV pulse profile. A comparison in absolute phase with the 
610 MHz radio profile shows alignment of the 7-ray pulses 
with two of three radio pulses. The luminosity of the pulsed 
emission (0.1-1 GeV) amounts L 7 = 1.64 • 10 34 • (AO/1 sr) • 
(d/5.7kpc) 2 erg s~ 4 which is ~ 7 % of the pulsar's total spin- 
down luminosity. The similarity of the X-ray and 7-ray pulse 
profile shapes of PSR J0218+4232 , and the apparent align- 
ment of the 7-ray pulses with two radio pulses at 610 MHz, 
bears resemblance to the well-known picture for the Crab pul- 
sar. This similarity, and the fact that PSR J0218+4232 is one 
of three millisecond pulsars (the others are PSR B 182 1-24 and 
PSR B 1937+21) which exhibit very hard, highly non-thermal, 
high-luminosity X-ray emission in narrow pulses led us to dis- 
cuss these millisecond pulsars as a class, noting that each of 
these has a magnetic field strength near the light cylinder com- 
parable to that for the Crab. None of the current models for 
7-ray emission from radio pulsars can explain the 7-ray spec- 
trum and luminosity of PSR J02 18+4232 . 

Key words: pulsars: individual: PSR J02 18+4232 - Stars: neu- 
tron - BL Lacertae objects: individual: 3C 66A - Gamma rays: 
observations - X-rays: stars 



1. Introduction 

Pulsed high-energy emission from millisecond (ms) pulsars has 
so far been detected at X-ray energies below ~ 10 keV for 



PSR B1937+21 (Takahashietal. 1999). The first two exhibit 



only five pulsars: P SR J0437-4715 ([Becker & Trumper 1993[ ), 
PSR J2124-3358 ( |Becker & Trumper 1999|), PSR B 182 1-24 
dSaitoetal. 1997|), PSR J0218+4232 flKuiper et al. 1998|) and 



broad X-ray pulses, have soft spectra and relatively low lu- 
minosities in the X-ray window, about 3 orders of magni- 
tude lower than derived for the latter three (L^~ 10 ~ 
10 32 erg s _1 assuming emission in a 1 sr beam). In addi- 
tion to the higher luminosity, these have very narrow X-ray 
pulses (intrinsically ~ 100/is or narrower) and hard power- 
law shape spectra measured up to ~ 10 k eV ( Saito et al. 1997 



Mineo et al. 200C , Takahashi et al. 1999 respectively), the two 
hardest spectra having indices as hard as ~ -0.65. This short 
observational summary suggests that this small sample can de 
devided in two distinct classes of ms pulsars: Class I, ms pul- 
sars with soft, low-luminosity X-ray emission in broad pulses; 
Class II, with highly non-thermal, high-luminosity X-ray emis- 
sion in narrow pulses. 

Millisecond pulsars not only differ from normal ra- 
dio pulsars in that their spin periods are 1 to 2 or- 
ders of magnitude shorter, reducing correspondingly their 
light cylinder radii, but particularly their surface mag- 
netic field strengths are 3 to 4 orders of magnitude 
weaker. Nevertheless, Bhattacharya & Srinivasan (1991) and 



[Stumer & Dermer (1994) showed that both of the competing 
classes of models for the production of 7-rays (polar cap and 
outer gap models) predict the production of detectable non- 
thermal emission up to the high-energy 7-rays for a sizable 
number of ms pulsars. An early systematic search for pulsed 
high-energy 7-ray emission from ms pulsars rendered, how- 



ever, only upper limits (Fierro 1995). In this paper we will 
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present circumstantial evidence for the first detection of pulsed 
high-energy gamma-ray emission from a Class II ms pulsar: 
PSR J02 18+4232. 



2 



Pulsed 7-ray emission from PSR J02 18+4232 



PSR J0218+4232 is a 2.3 ms pulsar in a two day or- 
bit around a low mass (~ 0.2 M Q ) white dwarf companion 



(Navarro et al. 1995; van Kerkwijk 1997). The dipolar perpen 



dicular magnetic field strength (B±) at the surface of the neu- 
tron star is estimated to be 4.3 x 10 8 G and the spin-down age 
is ~ 4.6 x 10 s years. The spin-down energy loss L sc j of the 
pulsar amounts ~ 2.5 x 10 35 erg s . The pulsar distance in- 
ferred from its dispersion measure and from the electron den- 



sity model of Taylor & Cordes (1993) is > 5.7 kpc 



Soft X-ray emission from the pulsar was first detected by 
Verbunt et al. (1996)|in a 20 ks ROSAT HRI observation. In a 



100 ks follow-up observation X-ray pulsations were discovered 
at a significance of about 5 a (Kuiper et al. 1998). The X-ray 
pulse profile is characterized by a sharp main pulse with an in- 
dication for a second peak at a phase separation of Ac/) ~ 0.47. 
The pulsed fraction inferred from the ROSAT HRI (0.1-2.4 
keV) data is 37 ± 13 %. It is interesting to note that also in the 
radio domain the source exhibits an unusually high unpulsed 
component of ~ 50 % ( [Mavarro et al. 1995 ). 

The ROSAT HRI provides no spectral information and the 
number of counts recorded in a far off-axis PSPC observation 
does not allow spectral modeling in the soft X-ray regime (0.1- 
2.4 keV). Also ASCA detected this source, however, the ob- 
servation was too short: no pulsation could be detected, and 
a spectral fit to the weak total excess resulted in a power-law 
photon index of -1.6 ± 0.6 (Kawai & Saito 1999). 

The spectral information for PSR J02 18+4232 improved 
enormously analyzing the data from a 83 ks BeppoSAX 
MECS (1.6-10 k eV) observation performed early 1999 
(Mineo et al. 2000). Pulsed emission was detected up to 10 



keV, the pulse profile clearly showing two peaks at the same 
phase separation of 0.47 which we reported earlier combining 
ROSAT HRI and PSPC observations ( [Kuiper et al. 1998[ ). The 
BeppoSAX MECS observation reveals that PSR J0218+4232 
exhibits the hardest pulsar X-ray spectra measured so far: Be- 
tween 1.6 and 10 keV one peak has a spectrum consistent with 
a power-law photon index of —0.84 and the other with an index 
of —0.42. The total pulsed spectrum can be described with an 

index -0.61 ( |Mineo et al. 2000) 

At high-energy 7-rays, Verbunt et al. (1996) noticed 
the positional coincidence of PSR J02 18+4232 with 
the second EGRET catalog source 2EG J0220+4228 
(Thompson et al. 1995), which was identified in the 
catalog and other publications with the BL Lac 



3C 66A ( 


Dingus et al. 1996 


Lamb & Macomb 1997 


). Using 



Mukherjee et al. 1997 



some additional EGRET 
observations, and applying a combination of spatial and 
timing analyses, Kuiper et al. (1999a) conclude that 2EG 
J0220+4228 is probably multiple: between 100 and 1000 
MeV PSR J02 18+4232 is the most likely counterpart, and 
above 1000 MeV 3C 66A is the best candidate counterpart. 
The third EGRET catalog ( Hartman et al. 199^ ), which is 
based on more viewing periods than the 2EG catalog, also 
identifies 3EG J0222+4253 (2EG J0220+4228) with 3C 
66A, rather than with the ms-pulsar. However, in a note on 
this source, they indicate that the identification with 3C 66A 



stems from the catalog position based on the > 1 GeV map. 
Furthermore, they confirm that for lower energies (100-300 
MeV) the EGRET map is consistent with all the source flux 
coming from the pulsar, 3C 66A being statistically excluded. 

In this paper we present the results of spatial, timing 
and pulse-phase resolved spatial analyses using all available 
EGRET (30 MeV - 30 GeV) data collected between Novem- 
ber 1991 and November 1998 in 5 observations with PSR 
J02 18+4232 within 25° of the pointing axis. Analysis of radio 
monitoring data of this pulsar provided us with an ephemeris 
valid over the total period of 7 years covering the EGRET 
observations, allowing phase folding of all selected EGRET 
events in a single trial. The resulting high-energy 7-ray pulse 
profile is compared with pulse profiles detected at X-ray ener- 
gies up to 10 keV, and in absolute phase with the radio profile 
at 610 MHz. The results are finaly discussed in relation to the 
Class II ms pulsars and the Crab, as well as with recent theoret- 
ical predictions for the production of X-ray and 7-ray emission 
in the magnetospheres of ms pulsars. 

2. Instrument description and observations 

EGRET (the Energetic Gamma Ray Experiment Telescope) 
aboard the Compton Gamma Ray Observatory (CGRO) has a 
(gas-filled) sparkchamber and is sensitive to gamma-rays with 
energies in the range 30 MeV to 30 GeV. In the mode used 
for most of the observations the field of view is approximately 
80° in diameter, although the instrument point-spread func- 
tion (PSF) and the effective area degrade considerably beyond 
30° off-axis. Its effective area is approximately 1500 cm 2 be- 
tween 200 and 1000 MeV, falling off at lower and higher ener- 
gies. The angular resolution is strongly energy dependent: the 
67% confinement angle at 35 MeV, 500 MeV and 3 GeV are 
10?9, 1?9 and 0?5 respectively. The energy resolution AE/E 
is ~ 20% (FWHM) over the central part of the energy range. 
Each registered event is time tagged by the on-board clock, 
serving also the other 3 CGRO instruments BATSE, OSSE 
and COMPTEL. The on-board time is converted to Coordi- 
nated Universal Time (UTC) with an absolute accuracy better 
than 100 /is, and a relative accuracy of 8 /is. For a continued 
proper sparkchamber performance regular gas replenishments 
of the sparkchamber are required in order to restore the effi- 
ciency after the gas has aged. The sparkchamber efficiency is 
therefore a function of time and energy. For a detailed overview 
of the EGRET detection principle and instrument characteris- 



tics, see Thompson et al. (1993). The inflight calibration and 



performance are presented in detail by Esposito et al. (1999). 
In this work we selected those CGRO Cycle I- VII Viewing 
Periods (VP) in which PSR J02 18+4232 , located at (Lb) = 
(139.508,-17.527), was less than 25° off-axis. In Table the 
details of each selected VP are given. 

COMPTEL, the imaging Compton Telescope aboard 
CGRO, is co-aligned with EGRET and had PSR J02 18+4232 
in its field of view during the same VP's as EGRET. COMP- 
TEL operates in the 0.75-30 MeV energy window and has an 
energy resolution of about 5-10 % FWHM, a large field of 



L.Kuiper et al. 

Table 1. EGRET observations used in this study with PSR J02 18+4232 less than 25° off-axis 
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VP# 


Start Date 


End Date 


Pointing 


direction 


Off-axis angle 


Eff.Exposure 


Sparkchamber efficiency 




TJD f 


TJD 


1(°) 


b(°) 


(°) 


(100-300 Me V; cm 2 s) 


(100-150 MeV)/(l-2 GeV) 


15.0 


8588.535 


8602.696 


152.75 


-13.40 


13.4 


3.209 x 10 8 


0.962/0.981 


211.0 


9043.646 


9055.631 


125.86 


-4.70 


18.5 


1.661 x 10 8 


0.870/0.935 


325.0 


9468.592 


9482.625 


147.08 


-9.06 


11.2 


2.512 x 10 8 


0.820 / 0.909 


427.0 


9951.603 


9967.581 


153.71 


-9.95 


15.7 


0.690 x 10 8 


0.269 / 0.632 


*728.7/9 


11078.646 


11120.603 


139.36 


-18.70 


1.2 


0.887 x 10 8 


0.180/0.250 



+ TJD = JD - 2440000.5 = MJD - 40000 

* EGRET in narrow field mode; opening angle FoV is 19° 



view (~ 1 steradian) and a location accuracy of ~ 1° (see 
Schonfelder et al. 1993). Events are time-tagged with a 0.125 



ms resolution. A timing analysis of PSR J02 18+4232 in the 
COMPTEL energy window did not yield a significant timing 
signal and subsequent imaging studies of the sky region con- 
taining PSR J0218+4232 did not show a source detection at 
the pulsar position. Therefore, only the flux upper limits are 
presented in this paper (see Sect. 8). 

OSSE, the Oriented Scintillation Spectrometer Experiment 
aboard CGRO, is a non-imaging detector system consisting of 4 
independent actively shielded NaI(Tl)-Cs(Na) phoswich detec- 
tors operating in the 0.05 to 10 MeV ener gy range, each havin g 
a 3? 8 x 11?4 (FWHM) field of view (see [lohnson et al. 1993| ). 
PSR J02 18+4232 was the primary target of OSSE during VP 
728.7 and VP 728.9 and was observed in event-by-event mode 
in the 50-150 keV energy band with a timing accuracy of 0. 125 
ms. Like in the case for COMPTEL, also OSSE did not detect 
a timing signature. Flux upper limits are given in Sect. 8. 



3. Spatial analysis 

Events arriving from within 30° off axis when EGRET is in 
full FoV mode and 19° in narrow field mode, are sorted in a 
3 dimensional data cube with galactic longitude, latitude and 
energy as axes. The longitude and latitude bin widths are 0?5, 
and 10 narrow "standard" energy ranges are selected: 30-50, 
50-70, 70-100, 100-150, 150-300, 300-500, 500-1000 MeV, 1- 
2, 2-4 and 4-10 GeV. 

Because the Earth atmosphere is the largest source of non- 
celestial 7-rays the events are subjected to an energy dependent 
zenith angle cut. We used the "standard" values for the 10 se- 
lected energy windows. The corresponding energy dependent 
exposure maps are calculated using the "exposure history" files 
taking into account the instrument calibration characteristics, 
the instantaneous timeline, the operation mode of the instru- 
ment and the time dependent spark chamber sensitivity factors 
(see ^sposito et al. 1999 and Table 1). 

To be consistent with the selection criteria used in the gen- 
eration of the exposure matrices we demand that the energy 
deposit in the TASC (Total Absorption Shower Counter) mea- 
sured by at least one of its PHA's is above a threshold of ~ 6.5 
MeV. 



The imaging method employed here is based on our Maxi- 
mum Likelihood Ratio (MLR) program, part of the CO MPTEL 
analysis software package COMPASS (de Vries 1994). In this 
program point sources are searched for on top of a diffuse back- 
ground model which describes the galactic and extra-galactic 
7-ray emission separately. The galactic component originating 
in cosmic -ray interactions with the protons of the atomic and 
molecular Hydrogen gas, as well as inverse Compton interac- 
tions of cosmic-ray electrons with the ambient photon field, is 
described by a combination of 2 different models: one results 
from the convolution of the EGRET PSF with the spatial distri- 
bution of the atomic Hydrogen column density and the second 
from the convolution with the spatial distribution of CO used 
as tracer for the molecular Hydrogen gas in the galaxy. The 
extra-galactic component is assumed to be isotropic. 

The image resulting from the Maximum Likelihood Ratio 
program is based on likelihood ratio tests at user defined grid 
points in a skyfield containing the object of interest. At each 
grid point (x s ky,Vsky) we determine the Maximum Likelihood 
under two hypotheses: 1) a description of the data in terms of 
the diffuse background models only (Ho) and 2) a description 
in terms of the diffuse background models and a point source 
at the (x s ky sky) position (Hi). Under the Hi hypothesis the 
number of counts (/i) expected in a measured sky pixel is 
given by : 



= <j-PSF 1 . 



.HI 



M? 1 + a co ■ + a Iso • M/, so (l) 
t j l j l j 



where M HI , M co and M Iso represent the convolved diffuse 
galactic and extra- galactic models. 

Because our mosaic of observations is composed of view- 
ing periods with pointing directions concentrated in a narrow 
band at low galactic latitudes (see Table [[]) the a Iso scale fac- 
tors are poorly constrained in the optimization process due to 
the dominating galactic components. So, we keep them fixed at 
values derived from a study of the extra-galactic 7-ray emis- 
sion using a much larger database containing all EGRET Cycle 



1,11 and III observations (Sreekumar et al. 1998 ) 



By optimizing the Likelihood under Hi with respect to 
its free scale parameters, er, a HI , a co we can derived the 
flux and flux uncertainty from a and its error for a putative 
source at position (x s ky,Usky)- From optimizations under Hi 
and Ho we can determine the Likelihood ratio A defined as 
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Fig. 1. MLR map for energies > 100 MeV of a sky re- 
gion centered on the EGRET source 2EG J0220+4228 / 3EG 
J0222+4253, combining data from the 5 observations listed in 
Table 0. The position of PSR J02 18+42 is indicated by a star 
symbol and of 3C 66A by a bullet. The contours start a 4ct de- 
tection significance level (1 d.o.f.) with steps of la. 



—2 ln(£-ft f) / C-Hx)- This quantity is distributed as a x 2 for 1 de- 
gree of freedom for a known source position and yields the 
source detection significance. 

The MLR map for energies > 100 MeV (Fig. [j]) con- 
firms the detection of the EGRET source 2EG J0220+4228 
/ 3 EG J0222+4253 (Thompson et al. 1995; Hartman et al. 
1999) at a <; 10er significance level for 1 degree of freedom, 
i.e. the source position is known. The Hi and 7io hypothe- 
ses include also the contributions from well-established 7-ray 
sources (Hartman et al. 1999) within a 30° radius around our 



target in order to describe the 7-ray sky near our target opti- 
mally. The binned event matrix for this integral energy win- 
dow is a combination of the matrices for the differential energy 
windows above 100 MeV, each with a different Earth zenith 
cut angle. The > 100 MeV exposure matrix is in this case a 
power-law weighted composition (index —2.1) of the differen- 
tial exposure matrices. This forms a consistent event/exposure 
set with respect to the applied selection criteria. 

We compared the derived optimum scale factors for 
the Galactic diffuse emission components with the find- 
ings from more detailed studies on this diffuse emission 



( |Strong & Mattox 1996j |Hunteret al. 1997| ) and found that our 
values are in all cases consistent with the published results. 

It is evident from Fig. [ljthat the high-energy 7-ray source is 
positionally consistent with both PSR J02 18+4232 and 3C 66 A 
(located at (l,b) = (140.143,-16.767)). The total excess con- 
tains 225 ± 27 counts. We analyzed this excess also in the dif- 
ferential energy windows: 100-300 MeV, 300-1000 MeV and 
1-10 GeV In each window the source was seen: 100-300 MeV 
^ 7. Oct detection significance and 138 ± 24 counts, 300-1000 
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-21.5 
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Fig. 2. MLR map showing 1, 2 and 3ct location confidence con- 
tours of the 7-ray source 2EG J0220+4228 / 3EG J0222+4253 
for 3 different energy windows. The shift of the excess to- 
wards the pulsar position for decreasing energies is evident. 
Between 100 and 300 MeV 3C 66A is located outside the So- 
contour, whereas between 1 and 10 GeV this is the case for 
PSR J02 18+4232. 

MeV ^ 7.0ct and 57 ± 12 counts and finally 1-10 GeV ^ 6.5ct 
and 22 ± 6 counts. The location confidence contours for the ex- 
cesses in the 3 broad energy windows are shown in Fig. [| This 
figure shows that 3C 66 A is the evident counterpart for the 1- 
10 GeV window (consistent with the third EGRET catalogue 
results ( [Hartman et al. 1999( ), whereas PSR J02 18+4232 is the 
most likely counterpart for the 100-300 MeV window. Between 
300 and 1000 MeV both sources contribute to the excess. 

For energies below 100 MeV we see indications for an 
excess, but the EGRET sensitivity is becoming too low and 
the spatial response too wide to draw firm conclusions. We 
estimated a 2cr flux upper limit for the spectrum of PSR 
J0218+4232 (see Sect. 8). 



4. Long-term time variability 

Earlier studies of 7-ray emission from spin-down powered pul- 
sars showed that they are steady 7-ray emitters (see e.g. the 



review by Thompson et al. 1997). On the contrary, most Ac 



tive Galactic Nu clei appeared to be hig hly variable at 7-ray 
energies (see e.g. Mukherjee et al. 1997 ). Therefore, we inves- 
tigated whether there is (absence of) evidence for time varabil- 
ity of 2EG J0220+4228 / 3EG J0222+4253, particularly for the 
1 00-300 MeV and 1 - 1 GeV bands, in which PSR J02 1 8+4232 
and 3C 66A appear to be the most likely counterparts, respec- 
tively. 

Using integration intervals of typically 2 or 3 weeks, the re- 
sults are shown in Fig. 0, The 100-300 MeV flux measurements 
are fully consistent with being constant, as expected for 7-ray 
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Table 2. Event extraction radius as a function of energy win- 
dow 



J.O 9000.0 9500.0 10000.0 10500.0 11000.0 
Time [TJD] 



VP 15 VP 211 VP 325 VP 427 



VP 728.7/9 _ 



9500.0 10000.0 
Time [TJD] 



Fig. 3. Long-term time variability of the 7-ray source 2EG 
J0220+4228 / 3EG J0222+4253 in different energy windows: 
100-300 MeV (top) and 1-10 GeV (bottom). The integration 
time for each data point is typically 2 or 3 weeks. The 100-300 
MeV flux points do not show time variability, while the 1-10 
GeV data points deviate at the ~ 2a level from being constant 
(variability index V = 1.33). Error bars are la; the shaded 
regions indicate the weighted mean ±la. 



emission from spin-down powered pulsars. The 1-10 GeV flux 
points show indications for variability and deviate at a ~ 2a 
level (93%) from being constant. According to the variability 
criteria defined by |McLaughlin et al. (1996) the 1 - 1 GeV vari- 
ability index V of 1.33 points to a variable nature of the 1-10 
GeV emission. This type of variability is indeed reminiscent 
of the behaviour observed frequently for the 7-ray emission 
from AGN. The above supports the conclusion from the spa- 
tial analysis, namely, that 2EG J0220+4228 / 3EG J0222+4253 
is multiple: above 1 GeV the BL Lac 3C 66A is the obvious 
counterpart, whereas below 300 MeV PSR J02 18+4232 is the 
most likely counterpart. 



5. Timing analysis 

In the timing analysis similar event selections have been ap- 
plied as in the spatial analysis, except we ignored the spe- 



Energy 


Extraction 


Enclosed 


window (MeV) 


radius (°) 


source fraction 


100- 150 


3.5 


0.53 


150- 300 


2.6 


0.56 


300- 500 


1.8 


0.56 


500-1000 


1.2 


0.56 



cific TASC (Thompson et al. (1993) I flags of the event triggers 
in the event selection process. Especially the selection on the 
TASC zero cross overflow bit (set to 1 if less than 6.5 MeV is 
deposited in the TASC), which is only effective for the lower 
energy 7-ray photons (< 150 MeV), is not taken into account. 
We verified this selection by a timing analysis of the Crab pul- 
sar (combining many Cycle — VI VP's) which showed that 
ignoring the TASC flags gives a significant improvement of the 
timing signal, particularly for energies below 100 MeV, with 
respect to the case in which we demand a TASC energy deposit 
of at least 6.5 MeV measured by one of its 2 PHA's. 

An additional difference in the selection procedure with the 
spatial analysis, where the spatial information of all events is 
used, is that we now have to specify an event extraction ra- 
dius around the pulsar position. Contrary to what is commonly 
used in the timing analysis of EGRET data, namely, select- 
ing events within an energy dependent extraction radius r ext 
of 5?85 • (E/100 MeV)~°- 53i containing approximately 67% 
of the counts from a point-source, with E the measured 7-ray 
energy, we optimized in each narrow energy window (e.g. 100- 
150 MeV) the signal-to-noise ratio S/N as a function of ex- 
traction radius taking into account the modelled (2d) spatial 
distribution of the optimized diffuse models and neighbouring 
sources as obtained in the spatial analysis (see e.g. the the- 



sis of Fierro 1995 p. 49-50). This method provides the optimal 
extraction radius for a given energy window and a given sky- 
background structure. The values obtained from this study for 
the narrow energy windows between 100 and 1000 MeV are 
listed in Table [| 

From our timing observations of PSR J02 18+4232 at radio 
wavelengths we obtained one single accurate ephemeris (rms 
error 85/is), which is listed in Table |[ The validity interval of 
this ephemeris covers almost 5 years and in view of the stable 
rotation behaviour observed for millisecond pulsars its validity 
should extend far beyond the indicated range. 

Phase folding the barycentered arrival times, taking into ac- 
count the binary nature of the system, of the selected events 
with measured energies between 100-1000 MeV from all ob- 
servations listed in Table [l] yields a 3.5cr modulation signif- 
icance applying a Z\ test ( [Buccheri et al. 1983 ) on the un- 
binned sample of pulse phases. An H-test ( de Jager et al. 198"^ ) 
in which the internal optimization of the number of harmon- 
ics is taken into account in the significance estimate yields a 
3.2cr modulation significance at an optimum number of har- 
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Table 3. Ephemeris of PSR J0218+4232 



Parameter 


Value 1 " 


Right Ascension (J2000) 


02 h 18 m 6" 350 


Declination (J2000) 


42° 32' 17'.' 44 


Epoch validity start/end (MJD) 


49092 - 50900 


Frequency 


430.4610674213 Hz 


Frequency derivative 


-1.4342 x 10" 14 Hzs" 1 


Epoch of the period (MJD) 


49996.000000023 


Orbital period 


175292.3020 s 


a ■ sin i 


1.98444 (lt-s) 


Eccentricity 





Longitude of periastron 





Time of ascending node (MJD) 


49996.637640 



* The last significant digit is given 




0.0 0.5 1.0 

Pulse phase 

Fig. 4. High-energy (100-1000 MeV) 7-ray pulse profile of 
PSR J02 18+4232 combining data from 5 separate viewing peri- 
ods. The modulation significance is ~ 3.5er applying a Z\ test. 
Background counts, e.g. from diffuse sky emission and possi- 
bly nearby sources, are included. The solid and broken lines 
indicate the Kernel Density Estimator (see text) with the ± la- 
uncertainty interval. A typical la error bar is shown. 



monies of 4. The 100-1000 MeV pulse profile is shown with 
10 bins in Fig, p] with superpo sed its Kernel Density Estima- 
tor (KDE; de Jager et al. 1986) with the ±lcr uncertainty inter- 



val. This KDE approaches the genuine underlying pulse profile 
(convolved with the instrumental time resolution) for an infi- 
nite number of events. The pulse profile shows one prominent 
narrow emission feature between phases ~ 0.6 and ~ 0.7 fol- 
lowing a broad less prominent pulse with maximum at phase 
~ 0.2. The phase separation of ~ 0.45 is remarkably similar to 
the value of ~ 0.47 observed at soft/medium energy X-rays by 
the ROSAT HRI (]Kuiperet al. 1998} ) and BeppoSAX MECS 
( Mineo et al. 2000| ; detailed comparisons will be presented be- 



low). 



We also produced phase distributions in broader differen- 
tial energy intervals. The pulse profiles for 100-300 MeV and 
300-1000 MeV both showed consistently the same narrow and 
broad pulses (Z\ probabilities 2.5er and 1.9<7, respectively). 
For 30-100 MeV and 1-10 GeV no hints for pulsation were 
found. 



6. Pulse phase resolved spatial analysis 

The pulse profile shown in Fig. ^ reaches a significance of 
~ 3.5cr, indicating that the probability is low, only 4.7 • 10~ 4 , 
that this deviation from a flat distribution is caused by a ran- 
dom fluctuation. Given the importance of the discovery of high- 
energy 7-ray emission from a millisecond pulsar, we investi- 
gated further whether there is additional support in our data to 
claim this detection. As explained above, for the timing analy- 
sis the events were selected within an extraction radius around 
the position of PSR J02 18+4232 using only - 56% of the 
events of a point source. In order to verify whether the source 
events outside the extraction radius (~ 44%) exhibit the same 
timing signature, we produced a pulse profile using all source 
events by performing a pulse phase resolved spatial analysis for 
energies between 100 and 1000 MeV. 

The procedure is the following: Construct a pulse profile 
by repeating the spatial analysis for events selected in different 
pulse phase intervals. Contrary to the phase folding we need 
to select the events in relatively broad phase intervals to have 
sufficient statistics to do the spatial analysis: We selected 10 
phase bins of width 0.1. 

In order to estimate first the contribution of 3C 66A, which 
is obviously independent of the pulsar phase, to the total high- 
energy 7-ray excess in the 100-1000 MeV energy band we have 
fitted this excess for the full [0,1] phase range in terms of point- 
sources at the positions of PSR J02 18+4232 and 3C 66 A. This 
yielded the following decomposition: the number of counts as- 
signed to PSR J0218+4232 and 3C 66A are 151 ± 52 and 
42 ± 51, respectively. The insignificant 3C 66A contribution, 
coming from events with energies > 300 MeV, is nevertheless 
taken into account as a small correction in the pulse phase re- 
solved spatial analysis (4.2 counts are assigned to 3C 66A for 
each 0.1 wide phase bin). Fitting then the measured 100-1000 
MeV spatial event matrices for each pulse phase slice in terms 
of a PSR J02 18+4232 model with a. free scale factor atop the 
galactic diffuse models (both with/ree scale factors), the (fixed) 
isotropic extragalactic component and all (fixed) nearby-source 
models including 3C 66A, we obtain the total number of counts 
correlating with a point-source at the PSR J02 1 8+4232 position 
for each phase slice. 

The resulting 10 bin pulse profile is shown in Fig. ||. The 
total number of source counts in this light curve is 153 (the 
background level = 0). Comparing Fig. || with the profile ob- 
tained from the timing analysis (Fig. Q), it is evident that the 
shape is statistically identical. For the phase folding we had 
selected only ~ 56% of the events for a real source (cf Table 
||). Scaling from the number of 153 source counts measured in 
Fig. |J to be consistent, the number of pulsar excess counts in 
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Pulsed 



Unpulsed 



Fig. 5. High-energy (100-1000 MeV) 7-ray pulse profile of 
PSR J02 18+4232 resulting from the pulse phase resolved spa- 
tial analysis. All background contributions are modelled out, 
including that of the nearby BL Lac 3C 66A. The profile is 
similar in shape to the profile from the timing analysis (Fig. 
The number of counts in this profile is a factor of ~ 1.8 higher 
than the excess counts in Fig. [|, as expected for a genuine pul- 
sar signal. A typical la error bar is given. 
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Fig. |] should be ~ 86, i.e. the backgound level should be at 
~ 22. It is evident from this comparison that the two profiles 
are fully consistent in shape as well as in number of counts 
in the timing signature. Thus, the pulsed signal is also present 
outside the dataspace confined by the used extraction radius, as 
expected for a real signal, i.e. the timing and spatial signatures 
are consistent with the detection of PSR J02 18+4232. 

A more well-known display of the same conclusion are 
"ON"-"OFF" maps, or "pulsed"-"unpulsed" maps. Guided by 
the shape of the 100-1000 MeV pulse profile in a 20 bin repre- 
sentation (see Fig. ^|e) we tentatively defined a "pulsed" phase 
interval as the combination of the phase ranges 0.05-0.40 and 
0.55-0.70 and an "unpulsed" interval as its complement. We 
then produced MLR maps selecting the events now also on 
their phase location in either of the 2 pulse phase windows for 
the 100-300, 300-1000 and 100- 1000 MeV energy ranges. The 
results are shown in Fig. |[ It is evident that the 100-300 MeV 
signal is confined within the "pulsed" interval, strengthening 
the conclusion that PSR J02 18+4232 is the counterpart of 2EG 
J0220+4228 for energies between 100 and 300 MeV. In the 
300-1000 MeV "unpulsed" MLR map ~ 4a residual emission 
is visible which can be explained by emission from 3C 66A 
and pulsed emission from PSR J02 18+4232 emitted outside 
the defined "pulsed" interval (e.g. possible contribution from 
a weak pulse near phase 0.9 in Figs. [|or |]). The overall picture 
for energies below 1000 MeV points to a very dominant PSR 
J02 18+4232 and a minor 3C 66 A contribution. 



Fig. 6. Pulse phase resolved MLR maps of the sky re- 
gion containing PSR J02 18+4232 in 3 different energy win- 
dows: 0.1-0.3 GeV, upper panels; 0.3-1 GeV, middle pan- 
els; 0.1-1 GeV, lower panels. Left, "pulsed" maps (phases: 
0.05 - 0.40 & 0.55 - 0.70); right, "unpulsed" maps (0.40 - 
0.55 & 0.70 — 1.05). The contours start at a 3a significance 
level in steps of la for 1 degree of freedom. PSR J02 18+4232 
is marked by a x and 3C 66A by a A symbol. The emission in 
the 100-300 MeV window is confined to the "pulsed" interval. 
For the 300-1000 MeV window the "unpulsed" interval shows 
~ 4cr residual emission. This can be explained by emission 
from 3C 66A in combination with pulsed emission from PSR 
J02 18+4232 not accounted for in the definition of the "pulsed" 
window. 

7. Multi-wavelengths profile comparisons 

7.1. Comparison with radio profiles 

The ephemeris of PSR J02 18+4232 given in Table [| and used 
for our 7-ray analysis, has been determined using Jodrell Bank 
observations at 610 MHz. The corresponding radio profile is 
shown in Fig. (see also Stairs et al. 1999). It is remarkable 
that the pulsar is practically never "off"; three pulses seem to 
cover the entire phase range from to 1 . 

Because the fiducial point in the 610 MHz radio profile 
defining the anchor point in the template used in the fitting 
process of the time of arrival of the radio pulses is known, its 
geocentric arrival time specified by the "Epoch of the period" 
in Table || can be translated to solar system barycentric arrival 
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time. This timestamp is subsequently converted to a phase zero 
taking into account the binary nature of the system. This phase 
zero value, corresponding to the fiducial point, is finally sub- 
tracted from the 7-event phases, obtained by the same fold- 
ing procedure, to align these with the radio profile. Thus, we 
can compare the 100-1000 MeV pulse profile in absolute phase 
with the 610 MHz radio profile. The aligned 7-ray pulse pro- 
file is shown in Fig. [7]e, now in 20 bins to allow a more detailed 
comparison. The bin width of ~ 115^s is comparable to the 
CGRO absolute timing accuracy of better than 100 /is. Also 
shown is the same KDE profile as shown in the 10 bin pulse 
profile in Fig. to aid the comparison of the two 7-ray his- 
tograms, given the low counting statistics. In order to guide the 
eye, the pulsar phases of the three maxima in the 610 MHz 
radio profile are indicated by vertical lines. 

I n Fig, f^b is also indicated the 1410 MHz radio pro- 
file (Kramer et al. 2000) which has been aligned by cross- 
correlation with the 610 MHz profile (phase uncertainty ~ 0.01 
in alignment). It is clear from this figure that the 2 emission 
features in the 7-ray pulse profile coincide within the absolute 
timing uncertainties with 2 of the 3 pulses in the 610 MHz ra- 
dio profile. Comparing the 610 and 1410 MHz radio profiles it 
is notable that one of these "radio/7-ray" pulses (at phase 0.62) 
coincides with a dip in the 1410 MHz profile, followed and pre- 
ceded by smaller pulses. Also between the two main emission 
features a shoulder is visible in the 1410 MHz profile which is 
absent in the 610 MHz one. 



7.2. Comparison with X-ray profiles 

We reported earlier significant detections of pulsed X-ray 
emission from PSR J02 18+4232 analysing ROSAT HRI data 
(4.8<7 modulation significance in the 0. 1-2.4 keV energy range; 



spectra of the two peaks (Mineo et al. 2000), make us confident 



Kuiperet al. 19980 and BeppoSA X MECS data (6.8a, 1.6-10 
keV energies; frlineo et al. 2000| ). In the BeppoSAX MECS 
analysis we used the same ephemeris of Table |] as in the 
present work. In the ROSAT HRI analysis, however, we used 
for the phase folding the extrapolated timing parameters from 
Navarro (1995). Given the availability of the new ephemeris 
which is valid over a nearly 5 year period and covers the 
ROSAT HRI observation, we decided for consistency reasons 
to reanalyze the 100 ks ROSAT HRI data. In addition, applica- 
tion of improved maximum likelihood algorithms in the spatial 
analysis to determine the centroid of emission in the X-ray map 
allowed for a better determination of the optimal extraction ra- 
dius (8")- The result is shown in Fig. ^c. The modulation sig- 
nificance has increased to 6cr (Z\ test), particularly the promi- 
nence of the second weaker pulse near phase 0.6 has improved 
in comparison with the result shown in Kuiper et al. (1998). 

The new ROSAT HRI profile can be compared with the 
BeppoSAX MECS profile (Fig. [7]d; |Mineo et al. 2000| ), which 
just overlaps in energy window. The alignment of the profiles 
was done by cross correlation, like in Mineo et al. (2000), since 
the uncertainties in the ROSAT and BeppoSAX absolute tim- 
ing are too large to allow an absolute comparison. The identical 
peak separations of ~ 0.47 and the consistent difference in the 



that the alignment is accurate. 

The next step is the alignment of the X-ray profiles with 
the absolute timing of the 7-ray and radio profiles. We cross 
correlated the most significant X-ray profile (from BeppoSAX 
MECS) with the EGRET profile, and applied the phase shift 
which corresponds to the highest probability in the correlation 
analysis to the aligned ROSAT HRI and BeppoSAX MECS 
profiles. These aligned profiles are shown in Fig. [7|. in which 
the BeppoSAX MECS and EGRET profiles are both displayed 
in 20 bins. It is obvious that all three high-energy profiles ex- 
hibit two pulses with the same phase separation of about 0.47. 
Fine structure in the gamma-ray profile, like the local maxi- 
mum at phase ~ 0.9, is not significant, even though the strong 
radio pulse at phase ~ 0.9 makes that phase "special". 

8. Multi-wavelength spectrum 

The X-ray spectrum of the pulsed emission from PSR 
J02 18+4232 between 1.6 and 10 keV is the hardest measured 
so far for any (millisecond) radio pulsar. The best power- 
law fit to the BeppoSAX MECS pulsed spectrum has an in- 
dex —0.61 ± 0.32. The spectrum becomes somewhat softer 
(index -0.94 ± 0.22) w hen a 27% DC component is in- 
cluded (Mineo et al. 2000). This DC component is visible in 



the ROSAT data flKuiper et al. 1998j ) and in the BeppoSAX 
data up to 4 keV. Above 4 keV there is no sign of a DC compo- 
nent. 

In the EGRET 7-ray data above 100 MeV, the signal seen 
from PSR J02 18+4232 is also consistent with being 100% 
pulsed. However, the detailed structure of the pulse profile is 
not clear, i.e., is there a phase interval in which the 7-ray signal 
is clearly off, or how wide are the wings of the pulses? Possi- 
bly, the pulsed 7-ray signal extends over the total phase range 
with only one or two very narrow dips, just like in the radio 
profile. Therefore, it is difficult to determine a background re- 
gion in the 7-ray pulse profile for the construction of a pulsed 
spectrum. We decided to determine the 7-ray spectrum using 
again the spatial maximum likelihood analysis, estimating the 
number of source counts (and then flux) on top of the diffuse 
background models and all relevant nearby sources, for the fol- 
lowing energy intervals: 30-100 MeV, 100-300 MeV, 300-1000 
MeV, 1-10 GeV The resulting flux values and upper limits are 
given in Table g for PSR J0218+4232 and the simultaneously 
derived values for 3C 66A in Table || (power-law photon in- 
dex ~ —1.5). Table [| also lists the upper limits derived for 
the simultaneous COMPTEL observations and the OSSE ob- 
servation during VP 728.7/9 (see Table |IJ). The COMPTEL 2a 
-upper limits are derived in a spatial analysis analoguous to the 
EGRET approach. The OSSE 2<r - upper limits are estimated 
from the statistically flat phase histograms according to the de- 
scription presented in Ulmer et al. (1991) assuming a duty cy- 
cle of 0.5. 

In Fig. U we have collected all available data for a total 
spectrum from radio up to high-energy 7-rays in the format 
E 2 x flux, showing the observed power per logarithmic energy 
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Table 4. CGRO flux estimates for PSR J0218+4232 



Energy 


Range 


Instrument 


Flux / 2cr upper limit 


[MeV] 




[ph / cm 2 s MeV] 


0.050 


0.073 


OSSE 




0.073 


0.103 




1 Q7 


0.103 


0.151 


OSSE 


< 0.87 ■ 10~ 3 


0.050 


0.151 


OSSE 


< 1.16 ■ 10~ 3 


0.75 


3 


COMPTEL 


< 2.87- 10" 5 


3 


10 


COMPTEL 


< 2.80 ■ 10~ 6 


10 


30 


COMPTEL 


< 2.76 ■ 10~ 7 


30 


100 


EGRET 


< 1.75 • 10 -8 


100 


300 


EGRET 


(7.71 ± 1.34) ■ 10~ 10 


300 


1000 


EGRET 


(3.86 ±2.08) ■ 10" 11 


1000 


10000 


EGRET 


< 1.25 ■ 10~ 12 



Table 5. CGRO EGRET time averaged flux estimates for 3C 
66A 



Energy Range 


Flux / 2a upper limit 


[MeV] 


[ph / cm 


sMeV] 


300 1000 


(4.24 ±2.11) 


• 10" 11 


1000 10000 


(3.38 ± 1.08) 


• io- 12 



interval. The very high luminosity at 7-ray energies between 
100 MeV and 1 GeV is striking and a large fraction of the total 
spin-down luminosity L s( j will be emitted in high-energy 7- 
rays. This fraction r/ ^ s can be estimated as follows: 



^obs — 



1.64 • 10 34 ■ {Afl/1 sr) • (d/5.7kpc)' 
2.36 • 10 35 • (J/10 45 gem 2 ) 



with AS1 the 7-ray beam size, d the distance to the pulsar and / 
the moment of inertia of the neutron star. Assuming Afi = 1 sr, 
d = 5.7 kpc and / = 10 45 gem 2 we obtain an efficiency of 
- 7% for PSR J02 18+4232. Over the 100-1000 MeV range the 
7-ray spectrum is soft and consistent with a photon power-law 
index of ~ —2.6. The extrapolation of the very hard spectrum 
between 0.1 and 10 keV is just in agreement with the OSSE 
upper limit(s). Fig. |8| suggests that the maximum luminosity is 
reached in the COMPTEL MeV range just below the COMP- 
TEL upper limits. 

9. Summary and discussion 

In this study we performed detailed spatial and timing analyses 
on PSR J0218±4232 using the high-energy 7-ray data from 
CGRO EGRET and found that we have good circumstantial 
evidence for the first detection of pulsed high-energy 7-rays 
from a Class II ms-pulsar, PSR J0218±4232 , namely: 

-1- The spatial distribution is consistent with the pulsar being 
detected: Between 100 and 300 MeV the EGRET source 



position is consistent with that of PSR J0218±4232 with 
the total signal concentrated in 2 pulses. The 100-300 MeV 
flux does not show time variability at a 2/3 weeks time 
scale, indicative for a steady 7-ray emitter like spin-down 
powered pulsars. Above 1 GeV the nearby (angular sepa- 
ration ~ 1°) BL Lac, 3C 66A, is the evident counterpart 
for the 7-ray excess. For energies between 300 MeV and 1 
GeV the pulsar and the BL Lac contribute to the excess. 
■2- Timing analysis (phase folding, using the timing parame- 
ters measured at radio wavelengths) in the 100-1000 MeV 
energy interval, selecting roughly 56% of the source pho- 
tons, yields a double-peaked pulse profile with a ~ 3.5<r 
modulation significance. The same pulsed signature is also 
present in the data outside the extraction radius used in the 
timing analysis, containing the remaining ~ 44% of the 
source photons. 

-3- The phase separation of ~ 0.45 of the two 7-ray pulses is 
similar to that measured between the two pulses at X-rays; 
a comparison in absolute time with the 610 MHz radio- 
profile shows alignment of the 7-ray pulses with two of the 
three radio pulses. 

EGRET detected six pulsars with overwhelming statisti- 
cal significance (Crab, Vela, Geminga, PSR B 1706-44, PSR 
B1951+32 and PSR B1055-52; see e.g. the review by Thomp- 
son et al. 1997). Compared to these six, the modulation signif- 
icance of PSR J0218±4232 falls only in the 3-4-a range, sim- 
ilar to the significance of the weak timing signals found with 
EGRET from P SR B0656+14 (|Ramanamurthy et al. 1996| ) and 
PSR B 1046-58 ( |Kaspi et al. 2000| ). The additional circumstan- 
tial evidence for the detection of PSR J0218±4232 , particu- 
larly the similarity of the double-peaked X-ray and 7-ray pulse 
profile shapes, and the fact that the X-ray spectrum measured 
for PSR J0218±4232 below 10 keV is the hardest measured 
for any pulsar (Mineo et al. 2000) increases the likelihood of 
the detection. Nevertheless, confirmation of the detections of 
PSR B0656+14, PSR B1046-58 and PSR J0218±4232 by fu- 
ture high-energy 7-ray missions like the Italian AGILE and 
NASA's GLAST is important. 

The nearby 3C 66A obviously complicated the analyses, 
but its contribution to the 7-ray excess in the skymaps has con- 
sistently been taken into account. The events detected from 
this BL Lac have no systematic effect on the double-peaked 
timing signature assigned to PSR J0218±4232 in the tim- 
ing analysis. However, our results show that earlier publi- 
cations on the spectrum of 3C 66A ( e.g. pingus et al. 1996 ; 



Mukherjee et al. 1997; Lin et al. 1999) should be revised, the 



time averaged spectrum is significantly harder than published 
earlier. 

InKuiperet al. (1998) and Mineo et al. (2000)the similarity 
of the double-peaked X-ray pulse profile of PSR J0218±4232 
with that of the Crab pulsar was noted and discussed. It is 
now striking that the observed 100-1000 MeV pulse profile 
of PSR J0218±4232 shows one narrow (~ 250/is) pulse pre- 
ceded ~ 0.45 in phase by a broader pulse, again a morphology 
very similar to that of the Crab pulsar 7-ray profile. The lat- 
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ter exhibits two distinct pulses at ~ 0.4 phase separation at 
X-ray and 7-ray energies, with the X-ray and 7-ray pulses be- 
ing aligned in absolute phase. Unfortunately, we cannot align 
the X-ray and 7-ray profiles of PSR J0218+4232 in absolute 
phase, but the similar phase separation suggests that the pulses 
are also aligned (see Fig. [7]). 

We noted in the Introduction that the surface magnetic field 
strengths of ms pulsars are 3 to 4 orders of magnitude weaker 
than that of normal radio pulsars. This makes the boundary 
condition for the production of 7-rays near the neutron star 
surface for ms pulsars much less favourable than for normal 
radio pulsars. It is, however, remarkable that the Crab pulsar 
and the members of the Class II ms pulsars have in common 
that the magnetic field strengths near the light cylinders B\ c 
are comparable (in the range (3 — 10) x 10 5 GauB). In fact, 
ranking all known radio pulsars by B\ c , the three Class II ms 
pulsars rank number 1 , 3 and 6, and Crab ranks number 2 (see 
also the discussion in Saito et al. 1997 ; Kuiper et al. 1998 and 
Takahashi et al. 19991). 



This is illustrated in Fig. ^, showing a scatter plot for all ra- 
dio pulsars of B\ c versus the spin-down flux, F s( j = E /(And 2 ), 
with E the total rotational energy loss rate and d the distance. 
The three Class II ms pulsars are clearly located at the ex- 
treme of the _Bi c distribution. The two Class I ms pulsars pos- 
sess significantly lower, more average values for ms pulsars. 
Also indicated are the 8 normal pulsars detected by EGRET in 
high-energy gamma-rays, as well as PSR B1509-58, dete cted 
by COMPTEL up to about 30 MeV ( |Kuiper et al. 1999b| ). As 
has been noted in earlier papers, F S( ± is a good indicator for the 
probability to detect hard X-ray and high-energy gamma-ray 
emission from normal radio pulsars. The only normal pulsar 
near the top of the F s( j distribution, not seen by EGRET is PSR 
B0540-69. This LMC pulsar is detected, however, at X-rays up 
to ~ 50 keV ( frjlmer et al. 1999 ). In order for ms pulsars to be 
seen with a hard X-ray spectrum (Class II), or even at high- 
energy gamma-rays (PSR J02 18+4232) a high value for B\ c 
seems to be required, in addition to a high F sc j. This suggests 
that _Bi c is a key parameter for models explaining the produc- 
tion of high-energy emission in the magnetospheres of ms pul- 
sars. 

Given in addition the similarities with the Crab of the 
high-energy pulse profiles (X-rays and also 7-rays for PSR 
J0218+4232 ) this suggests that the pulsed high-energy non- 
thermal emission from the Class II ms pulsars and the Crab 
pulsar have a similar origin in the pulsar magnetosphere, quite 
likely in a vacuum gap near the light cylinder. We know from 
radio observations, however, that the Crab has an orthogonal 
alignment, while PSR J0218+4232 is a nearly aligned rotator 
(Navarro et al. 1995, Stairs et al. 1999). Unfortunately, a pa- 
rameter which is also important in this discussion on the geom- 
etry, the impact angle, has only been determined with large un- 
certainties, and therefore the line-of-sight information for PSR 
J0218+4232 is unconstrained (Stairs et al. 1999). 

If indeed, X-ray emission and 7-ray emission from Class II 
ms pulsars has to be produced in a vacuum gap near the light 
cylinder, the vacuum gap has to be very short in order to have 



narrow and aligned pulses at X-rays and 7-rays, given the very 
strong curvature of the magnetic field lines in ms-pulsar mag- 
netospheres. In addition, the potential drop has to be very high 
over this short length to accelerate the particles to the ener- 
gies required for high-energy gamma-ray production. It is ob- 
vious that continuous acceleration of particles and production 
of X-rays and 7-rays from the surface of the neutron star along 
the curved magnetic field lines till the light cylinder radius (for 
PSR J0218+4232 only 111 km) will not render the narrow and 
aligned pulses at X-rays and 7-rays. 

The Crab pulsar has also its two X-ray and 7-ray pulses 
aligned in absolute phase with two of the three radio pulses, 
leading to a consistent picture in which the high-energy pulses 
and the aligned radio pulses are produced in the same zones in 
the magnetosphere (see e.g. feomani & Yadigaroglu 1995 ). The 
apparent alignment of the 7-ray pulses of PSR J0218+4232 
with two of three pulses measured at 610 MHz suggests also 
that some of the radio pulses are produced in the same zones in 
the magnetosphere as the 7-ray pulses. However, we would first 
like to see a better radio estimate of the viewing angle for the 
PSR J0218+4232 system, and a confirmation of the absolute 
alignment using new and better observations at X-ray energies, 
before making further speculations on the geometry. 

Theoretical models attempting to explain the high-energy 
electro-magnetic radiation from spin-down powered pulsars 
are divided in two main catagories distinguished by the pro- 
duction sites of the radiation in the pulsar magnetospheres. 
The first class of models, polar cap (PC) models, rely on 
the acceleration of charged particles along the open field 
lines near the magnetic pole(s) followed by cascade processes 
given rise to high-energy electro-magnetic radiation (see e.g 
paugherty & Harding 1994 , 1996). In the second class of mod- 
els, outergap models (OG), the acceleration of charged parti- 
cles and subsequent generation of high-energy radiation takes 
place in vacuum gaps near the pulsar light cylinder (see e.g. 



Cheng et al. 1986a,b and Ho 1989). Unfortunately, for the case 



of ms pulsars no detailed self-consistent model calculations ex- 
ist for either class of models, allowing predictions for different 
observational aspects, e.g. pulse phase resolved spectra, pulse 
shapes, efficiencies. In most cases only one aspect of the emer- 
gent high-energy radiation is addressed. 

The PC model elaborated by a Polish group (Bulik, Dyks 
& Rudak), for example, only focusses on the emergent high- 
energy electro-magnetic spectrum from ms pulsars from X-rays 
up to high-energy 7-rays, while the pulse shape is ignored. This 
group predicts a dominating Synchrotron component over the 
entire X-ray /soft 7-ray band (0. 1 ke V - 1 MeV) with a spec- 
tral photon index of —1.5 ( pyks & Rudak 19"99 ). This does 
not agree with the much harder photon indices of il — 1 ob- 
served for PSR J02 18+4232, PSR B 1937+21 and PSR B 1821- 
24. The predicted 7-ray spectrum, dominated by curvature ra- 
diation, peaks between 10 GeV and 100 GeV and even an in- 
verse Compton scattering component is predicted at TeV ener- 



gies (Bulik & Rudak 1999; Bulik et al. 2000). The maximum 



in the observed spectrum of PSR J0218+4232 (vF v or E 2 F 
representation) is located in the 1-100 MeV range, also in con- 
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tradiction with their model prediction (see Fig. ||). Their 7-ray 
flux prediction for PSR J02 18+4232 is even more than a factor 
of 10 below the expected GLAST sensitivity level, thus not at 
all detectable by the less sensitive EGRET telescope for which 
we present the results. 



The polar cap cascade model of Zhang & Harding (2000) 



including now also, compared to earlier versions, inverse 
Compton scattering of higher generation cascade pairs pro- 
vides predictions for both the X-ray and 7-ray luminosities 
of spin-down powered pulsars, including ms pulsars. In the 
soft/medium energy X-ray band the model predicts a thermal 
origin of the spectral features of the pulsed emission from 
ms pulsars. This is inconsistent with the observed non-thermal 
(very) hard pulsed spectra of the 3 Class II ms pulsars. How- 
ever, for the Class I ms pulsars this could be in agreement with 
the observed spectral properties. 

Zhang and Harding also predict that ms pulsars usually 
have a considerable high-energy 7-ray luminosity, but due to 
their weak magnetic field strengths, resulting in quite high pho- 
ton escape energies, the emergent 7-spectrum is very hard. The 
latter is not in agreement with the observed soft high-energy 
(photon Power-law index of ~ —2.6 for energies between 100 
MeV and 1 GeV) 7-ray spectrum of PSR J02 18+4232 . 

Thus, so far the PC scenario based models appear to be 
unsuccessful in explaining the observed X-ray and 7-ray prop- 
erties of the Class II ms pulsars. 

An outergap model aiming at predicting pulsed and un- 
pulsed 7-ray em ission from ms pulsars was presented by 
Wei et al. (1996) . This model predicts a spectral photon index 
of —2 for the pulsed emission from a ms pulsar for the en- 
ergy range of ~ 10 keV to ~ 500 MeV, in contradiction with 
the spectrum we show in Fig. | for PSR J02 18+4232 . The 
model predicts also a harder unpulsed component with a spec- 
tral photon index of —1.5, dominating the pulsed component 
above ~ 500 MeV. We have not detected this component for 
PSR J0218+4232 at energies above 100 MeV. 

Concerning the energetics of the 7-ray emission of PSR 
J0218+4232 it is interesting to compare the observed 7-ray 
efficiency (fraction of the total spin-down luminosity) 
of ~ 0.07 with theoretically derive d efficiencies. For the 
PC model of Zhang & Harding (2000) the efficiency scales as 
rjpQ oc p ■ t - 5 with p the pulse period and r the characteris- 
tic age of the pulsar. Expressed in the Crab pulsar efficiency 
VCmb we find for PSR J0218+4232 that 7? 0218 - 45 x 7y Crab , 
which translates to an efficiency of ?7o218 ~ 0.05 substituting 
the measured Crab 7-ray efficiency of about 0.001. The thick 
OG model of Zhang & Cheng (1998^ yields the following ex- 
pression for the 7-ray efficiency: t/qq oc p 2 ■ t 6 / 7 . This trans- 
lates to ?7o218 ~ 300 x rjQ la fo, which means that t?0218 ~ 0.33. 
Thus, within the framework of both PC and OG models the ex- 
pected 7-ray conversion efficiency is very high, approximately 
in accordance with the measured efficiency of about 0.1. How- 
ever, it should be noted that both models predict an even higher 
efficiency for e.g. PSR J0437-4715, a Class I ms pulsar. This 
pulsar is very nearby but has not been detected as a 7-ray 
source/pulsar ( Fierro et al. 1995 ). 



The circumstantial evidence presented in this paper for the 
detection of pulsed high-energy 7-rays from ms pulsar PSR 
J02 18+4232 opens a new window in the study of the magneto- 
spheric properties of spin-down powered pulsars. It is unfortu- 
nate that we cannot repeat this observation with EGRET any- 
more. Therefore, deep searches for high-energy 7-ray emission 
from the Class II ms pulsars with future more sensitive gamma- 
ray missions like GLAST and AGILE are very important. But 
also earlier sensitive observations at the harder X-rays above 
10 keV are very important to bridge the observational gap. Par- 
ticularly the ES A mission INTEGRAL might be able to extend 
the hard spectra measured below 10 keV to as high as a few 
MeV. 
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Fig. 7. Multi-wavelength pulse profiles of PSR J02 18+4232. 
Radio pulse profile at 610 MHz and 1410 MHz are shown in 
panels a and b respectively, c) reanalyzed ROSAT HRI 0.1- 
2.4 keV pulse profile, d) BeppoSAX MECS 1.6-10 keV pulse 
profile and e) 100-1000 MeV EGRET pulse profile. Indicated 
as dotted lines are the positions of the 3 pulses in the 610 MHz 
radio profile. The X-ray profiles are aligned at their highest cor- 
relation value with the EGRET (absolute) 100-1000 MeV pulse 
profile. Typical ±ler error bars are indicated in the X and 7-ray 
profiles. 
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Fig. 8. Multi-wavelength spectrum of the pulsed emission from 
PSR J0218+4232. The high-energy spectrum is characterized 
by a rapid rise at X-rays, followed by a flattening in the MeV- 
regime and a decline at high-energy 7-rays. Maximum lumi- 
nosity is reached in the MeV domain, however, the peak flux is 
just below COMPTEL's current sensitivity level. The EGRET, 
COMPTEL and OSSE upper limits are 2a. Error bars: ±la. 
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Fig. 9. Magnetic field strength near the light cylinder B\ c ver- 
sus pulsar spin-down flux F S( j for the sample of normal radio 
pulsars (□) and ms pulsars (o). The eight (see text) normal pul- 
sars detected above 1 MeV are indicated by a filled square sym- 
bol; Class I ms pulsars by an encircled filled circle; Class II by 
a filled circle. 



